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INTRODUCTION: 


The  majority  of  patients  with  advanced  stage  epithelial  ovarian  cancer  (EOC)  present  with 
advanced  stage  disease,  which  is  currently  treated  by  cytoreductive  surgery  and  chemotherapy. 
Approximately  10%  of  EOC  patients  cannot  be  successfully  cytoreduced  by  surgery  and  20%  are 
intrinsically  resistant  to  chemotherapy  or  develop  chemoresistant  disease  within  one  year  from 
initial  treatment.  Currently,  ovarian  cancer  surveillance  and  subsequent  therapies  are 
implemented  on  a  “watch-and-wait”  basis  because  there  are  no  reliable  biomarkers  to  identify 
patients  with  adverse  outcomes  on  standard  treatment.  To  identify  biomarkers  that  predict 
adverse  outcome  in  patients,  we  studied  the  key  processes  involved  in  metastatic  ovarian  cancer 
progression,  including  changes  in  the  tumor  microenvironment.  This  led  to  the  identification  of  a 
stromal/extracellular  matrix  gene  signature  that  correlates  with  poor  patient  survival.  In  the  first 
part  of  this  project,  we  have  identified  and  optimized  gene  signatures  for  the  identification  of 
patients  who  are  unlikely  to  benefit  from  standard  surgery  and/or  chemotherapy.  Since  most  of 
the  signature  genes  were  expressed  in  the  cancer  stroma  (fibroblasts  and  extracellular  matrix 
secreted  by  fibroblasts),  we  focused  on  the  molecular  characterization  of  cancer  fibroblasts. 
Cancer  fibroblasts  actively  contribute  to  ovarian  cancer  progression  and  are  viewed  as  a 
promising  therapeutic  target.  However,  the  role  of  cancer  fibroblasts  as  collaborators  in  cancer 
progression  has  recently  been  challenged.  Ablation  of  a-smooth  muscle  actin  (aSMA)-positive 
fibroblasts  in  a  mouse  model  of  pancreatic  cancer  has  resulted  in  increased  cancer  progression 
and  metastasis,  leading  authors  to  conclude  that  fibroblasts  have  a  preventive  role  in  the 
progression  of  some  cancers  [1].  Another  recent  study  in  a  mouse  model  of  pancreatic  cancer 
concluded  that  cancer  genotype  determines  the  activation  state  of  cancer  fibroblasts  [2]. 

The  contradictory  results  in  different  cancer  models  could  be  explained  by  different  roles  of 
fibroblasts  in  different  cancer  types,  i.e.  fibroblasts  could  be  promoting  ovarian  cancer  and 
inhibiting  pancreatic  cancer.  Alternatively,  in  all  cancer  types  fibroblasts  prevent  cancer 
progression  until  they  receive  activating  signals  from  cancer  cells  and  convert  into  ‘activated 
fibroblasts’,  which  in  turn  confer  invasive  and  metastatic  abilities  upon  cancer  cells  [3], 
Therapies  that  target  all  fibroblasts  are  counterproductive  and  likely  to  result  in  the  death  of 
normal  fibroblasts  and  significant  toxicity. 

Preferential  targeting  of  activated  fibroblasts  has  been  challenging  because  activated  fibroblasts 
are  poorly  understood  at  the  molecular  level.  During  activation,  fibroblasts  exhibit  phenotypic 
changes  that  partially  overlap  with  myofibroblastic  changes  during  wound  healing, 
inflammation,  and  fibrosis,  including  secretion  of  specific  ECM  components,  cytokines  and 
growth  factors  [4,  5].  Several  markers  have  been  used  to  distinguish  activated  from  non-activated 
fibroblasts:  a-smooth  muscle  actin  (aSMA,  encoded  by  gene  ACTA2),  fibroblast  activation 
protein  (FAP),  podoplanin  (PDPN),  palladin  (PALLD),  tenascin-C  (TNC),  platelet-derived 
growth  factor  receptor  a  (PDGFRa),  and  NG2  chondroitin  sulfate  proteoglycan  (CSPG4). 
However,  these  markers  are  frequently  expressed  in  other  cells  within  the  cancer  stroma,  such  as 
vascular  smooth  muscle  cells,  pericytes,  and  mesenchymal  stem  cells.  This  lack  of  specificity 
could  pose  problems  in  therapeutic  targeting  and  underscores  the  need  to  better  understand  the 
molecular  characteristics  of  activated  fibroblasts  in  order  to  develop  more  precise  and  less  toxic 
targeted  therapies. 
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We  previously  identified  COL11A1  as  part  of  gene  signatures  associated  with  suboptimal 
debulking,  rapid  recurrence  and  poor  overall  survival  in  ovarian  cancer  [6,  7],  In  situ 
hybridization  and  immunohistochemistry  in  ovarian  cancer  revealed  that  COL11A1  mRNA  and 
pro-protein  are  primarily  expressed  in  cancer  fibroblasts  [7,  8],  We  explored  the  suitability  of 
COL11A1  as  a  pan-cancer  marker  of  activated  fibroblasts  and  used  it  as  a  ‘seed’  to  identify  the 
transcription  signature  of  activated  fibroblasts  in  13  epithelial  cancer  types.  We  show  that  the 
COL11A1- coexpressed  gene  set  is  highly  conserved  in  these  13  cancer  types,  indicating  that  the 
fibroblast  reaction  to  cancer  cells  is  independent  of  the  organ  site-of-origin  and  of  the 
transforming  events  within  cancer  cells. 


2.  KEYWORDS: 

Ovarian  cancer,  tumor  microenvironment,  gene  signature,  carcinoma-associated  fibroblasts, 
clinical  outcome 

3.  ACCOMPLISHMENTS: 

■  What  were  the  major  goals  of  the  project? 


Specific  Aim  1  (specified  in  proposal) 

Timeline 

%  Completed 

Major  Task  1 

Identify  gene  signatures  for  the  prediction  of  poor  outcome 

Months 

Subtask  1 

Select  biologically  relevant  covariates  and  build  a  multivariate  model  for  the  analysis 
of  3  datasets  (TCGA,  n=403;  GSE26712,  n=185;  and  GSE51088,  n=122;  these  are 
public  datasets  with  de-identified  patient  information). 

1-3 

100% 

Subtask  2 

Analyze  datasets  individually  and  in  combination;  derive  gene  signatures  using 
multiple  statistical  methods;  correlate  with  overall  survival,  progression-free 
survival,  residual  disease  and  other  outcomes 

4-6 

100% 

Subtask  3 

Identify  small  subsets  of  predictive  genes  and  their  interactions 

7-8 

100% 

Milestone  Achieved 

A  gene  signature  consisting  of  8-15  genes  with  high  predictive  power  in  all  three 
datasets 

8 

100% 

Major  Task  2 

Optimize  the  gene  signature  for  the  prediction  of  poor  outcome 

Subtask  1 

Assess  the  predictive  accuracy  of  the  gene  signatures  using  independent  datasets 
(GSE9891  and  GSE3149) 

9-10 

100% 
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Subtask  2 

Generate  a  test  qPCR  set  from  frozen  samples  of  patients  with  extreme  outcomes  (10 
with  <1  year  survival  and  10  with  >7  year  survival);  validate  by  qPCR  up  to  30 
prioritized  genes  that  have  extremely  high  predictive  power  in  individual  datasets  but 
are  not  present  in  all  3  datasets 

11-12 

10% 

Subtask  3 

Validate  gene  signature  accuracy  using  statistical  methods 

13-14 

0% 

Milestone  Achieved 

An  independently  validated  set  of  8-15  genes  with  high  predictive  power 

14 

0% 

Major  Task  3 

Validate  the  gene  signature  for  the  prediction  of  poor  outcome 

Subtask  1 

Identify  and  collect  200  primary  ovarian  cancer  patient  samples  with  annotated 
demographic,  pathologic  and  clinical  information  and  follow-up  (all  patient  samples 
will  be  de-identified) 

3-7 

100% 

Subtask  2 

Cut  and  stain  slides  (1  H&E  +  9  unstained  sections);  evaluate  the  suitability  of  each 
sample  by  pathologic  examination  of  tumor  sections  and  circle  the  area  on  the  slide 
for  RNA  isolation 

8-9 

100% 

Subtask  3 

Isolate  RNA  from  slides;  perform  quality  control 

9-12 

100% 

Subtask  4 

Design  the  Nanostring  assay  for  the  quantification  of  the  signature  genes;  perform 
quality  control;  collect  data 

13-15 

100% 

Subtask  5 

Analyze  data  using  statistical  methods;  correlate  with  overall  survival,  progression- 
free  survival,  residual  disease  status  and  other  outcomes  that  could  be  used  to 
improve  clinical  management  of  ovarian  cancer  patients 

16-21 

0% 

Subtask  6 

Validate  the  gene  signature  assay  using  statistical  models  and  risk  prediction  models 
with  known  parameters 

21-22 

0% 

Subtask  7 

Submit  a  manuscript  on  the  predictive  power  of  the  optimized  gene  signature  using 
microarray  data  and  Nanostring  assay  data  and  deposit  RNA  expression  data  into 
public  repository  (GEO) 

Plan  an  academic  multi-center  validation  of  the  Nanostring  signature  assay  as 
required  prior  to  FDA  validation 

23-24 

0% 

Milestone  Achieved 

Validated  gene  signature  gene  assay  for  the  prediction  of  clinical  outcome(s)  using 
paraffin-embedded  tumor  tissues 

22 

0% 
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What  was  accomplished  under  these  goals? 

1 )  major  activities 

In  the  first  part  of  this  project,  we  have  identified  gene  signatures  for  the  identification  of 
patients  who  are  unlikely  to  benefit  from  standard  surgery  and/or  chemotherapy.  Since  most  of 
the  signature  genes  were  expressed  in  the  cancer  stroma  (fibroblasts  and  extracellular  matrix 
secreted  by  fibroblasts),  we  focused  on  the  molecular  characterization  of  cancer  fibroblasts. 

2 )  specific  objectives 

To  effectively  target  fibroblasts  in  ovarian  cancer,  it  is  necessary  to  accurately  characterize  the 
molecular  characteristics  of  the  activated  fibroblast  populations  in  different  cancer  types  and 
determine  whether  ovarian  cancer  activated  fibroblasts  are  different  from  activated  fibroblasts  in 
other  cancer  types.  Considering  the  artificial  nature  and  limitations  of  mouse  cancer  models,  it  is 
important  to  conduct  this  analysis  in  human  tissues. 

Our  objective  was  to  conduct  an  in-depth  analysis  of  public  expression  datasets.  First,  we 
defined  COL11A1  as  a  highly  specific  biomarker  of  activated  fibroblasts  in  multiple 
epithelial  cancer  types.  Second,  using  COL11A1  as  a  ‘seed’  to  identify  coexpressed  genes, 
we  demonstrated  that  activated  fibroblasts  express  a  highly  conserved  gene  signature  across 
genetically  different  epithelial  cancer  types. 

3)  significant  results  or  key  outcomes,  including  major  findings,  developments,  or  conclusions 

3.3.1.  COL11A1  is  associated  with  cancer  progression  and  adverse  clinical  outcomes  in 
ovarian  cancer 

COL11A1  mRNA  expression  has  been  associated  with  poor  survival  in  ovarian  cancer  [7,  9].  To 
elucidate  the  underlying  biology  that  could  result  in  poor  survival,  we  investigated  its  expression 
in  ovarian  and  colon  cancers.  Using  a  comprehensively  annotated  microarray  database  for  3431 
human  ovarian  cancers  [10],  we  show  that  increased  expression  of  COL11A1  mRNA  is 
associated  with  disease-specific  and  disease-free  survival  (Fig.  1A)  as  well  as  with  clinical  and 
molecular  parameters  such  as  increased  cancer  stage  and  grade  and  mesenchymal  molecular 
subtype  (Fig.  IB). 
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Disease-Specific  Survival 


HR  =  0.761 3  (0.6720  -  0.8624) 
Log-rankp  =  0 

COL1 1 A1  low  high 

Median  Survival  52.77  42.17 
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Log-rankp  =  0 
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Disease-Free  Survival 
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Median  Survival  25.00  18.10 


Disease-Free  Survival 
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Fig.  1.  COL11A1  expression  is  associated  with  adverse  clinical  parameters.  (A)  Kaplan- 
Meier  survival  plots  and  (B)  plots  of  COL11A1  expression  in  individual  subsets  of  ovarian 
cancer  were  generated  using  the  ovarian  microarray  gene  expression  database  CSIOVDB 
(csibio.nus.edu.sg/CSIOVDB/CSIOVDB.html). 

3.3.2.  COL11A1  is  expressed  in  a  subset  of  aSMA-positive  cancer  fibroblasts  and  can  be 
induced  in  normal  fibroblasts  by  the  presence  of  cancer  cells 

To  determine  if  COL11A1  expression  is  associated  with  fibroblast  activation,  we  used  aSMA  as 
a  marker  of  activated  fibroblasts  [11],  Comparison  of  aSMA  immunohistochemistry  and 
COL11A1  in  situ  hybridization  in  a  tissue  microarray  consisting  of  primary,  metastatic  and 
recurrent  ovarian  cancers  from  42  patients  showed  that  COL11A1  is  expressed  in  a  subset  of 
aSMA-i-  fibroblasts  (Fig.  2A).  Unlike  aSMA,  COL11A1  was  not  expressed  in  blood  vessels  (red 
arrows)  or  in  fibroblasts  surrounding  the  cancer  (blue  arrows)  (Fig.  2A).  In  sections  of  metastatic 
ovarian  cancer,  we  observed  that  COL11A1  -positive  cells  are  confined  to  the  intratumoral  and 
immediate  peritumoral  fibroblasts  (Fig.  2B),  suggesting  that  COL11A1  expression  may  be 
induced  by  cues  received  from  epithelial  cancer  cells.  To  test  if  cancer  cells  can  induce 
COL11A1  expression  in  fibroblasts,  we  co-cultured  immortalized  normal  ovarian  fibroblasts 
(INOFs)  with  three  different  ovarian  cancer  cell  lines  (OVSAHO,  OVCAR3,  and 
KURAMOCHI).  COL11A1  expression  in  INOFs  was  most  strongly  induced  by  direct  co-culture 
with  ovarian  cancer  cell  lines  although  weak  induction  occurred  by  indirect  co-culture  on  a 
Transwell  membrane  (Fig.  2C).  The  induction  of  COL11A1  in  fibroblasts  in  the  presence  of 
cancer  cells  was  confirmed  by  analysis  of  the  public  expression  dataset  GSE52104  in  which  two 
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types  of  presumptive  cancer-associated  fibroblast  precursor  cells,  mesenchymal  stem  cells 
(MSCs)  and  immortalized  normal  ovarian  fibroblasts  (INOFs),  were  either  cultured  alone  or  co¬ 
cultured  with  normal  ovarian  surface  epithelial  cells  (IOSE)  or  epithelial  ovarian  cancer  cells 
(EOC)  using  a  Transwell  membrane  [12].  COL11A1  mRNA  was  statistically  significantly 
upregulated  when  MSCs  and  INOFs  were  co-cultured  with  EOC  but  not  IOSE  (Fig.  2D), 
indicating  that  cancer  cells  have  a  greater  capacity  than  normal  cells  to  induce  COL11A1 
expression  in  fibroblasts. 
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Fig.  2.  COL11A1  is  expressed  in  cancer-activated  fibroblasts.  (A)  Comparison  of  aSMA 
immunohistochemistry  and  COL11A1  in  situ  hybridization  in  a  metastatic  ovarian  cancer  sample. 
Red  arrows  indicate  blood  vessels.  Blue  arrows  indicate  fibroblasts  surrounding  the  tumor 
nodule.  (B)  Distribution  of  COL11A1- positive  fibroblasts  in  relation  to  cancer  cells.  The 
estimated  percent  of  COL11A1- positive  fibroblasts  is  shown  on  the  right.  The  image  is 
representative  of  metastatic  and  recurrent  ovarian  cancer  samples,  which  typically  express  higher 
levels  of  COL11A1  than  primary  ovarian  cancers.  (C)  Quantitative  RT-PCR  levels  of  COL11A1 
in  sorted  (FACS)  immortalized  normal  ovarian  fibroblasts  (INOFs)  grown  alone  or  co-cultured 
with  ovarian  cancer  cell  lines  (OVSAHO,  OVCAR3,  KURAMOCHI)  that  were  either  separated 
from  INOFs  by  a  Transwell  membrane  or  directly  mixed  with  INOFs.  Statistical  analyses  were 
performed  between  INOFs  grown  alone  and  INOFs  co-cultured  with  ovarian  cancer  cells 
(*p<0.05;  **p<0.01;  ***p<0.001;  ns,  not  significant).  Error  bars  indicate  standard  deviation.  (D) 
Levels  of  COL11A1  in  the  GSE52104  expression  dataset  in  which  mesenchymal  stem  cells 
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(MSCs)  or  immortalized  normal  ovarian  fibroblasts  (INOFs)  were  either  cultured  alone  or  co¬ 
cultured  with  IOSE4  normal  epithelial  cells  (IOSE)  or  HEYA8  epithelial  ovarian  cancer  cells 
(EOC)  using  a  Transwell  membrane.  Inverse-log2  data  of  different  COL11A1  probes  were 
averaged,  then  log2-transformed.  The  data  were  extracted  for  statistical  analyses  using  GraphPad 
Prism  6.  Data  are  represented  as  the  mean  ±  SEM.  Intergroup  differences  were  assessed  by  the 
Student’s  t-test.  *p<0.05;  ****p<0.0001. 


3.3.3.  COL11A1  is  among  the  most  differentially  expressed  genes  between  cancers  and 
corresponding  benign  tissues 


Comparison  of  COL11A1  expression  in  17931  cancers  and  3503  normal  tissues  (U133Plus2 
platform)  and  9258  cancers  and  4087  normal  tissues  (U133A  platform)  using  the  Gene 
Expression  across  Normal  and  Tumor  tissue  (GENT)  portal  [13]  revealed  that  COL11A1  mRNA 
is  elevated  in  most  cancers  in  comparison  to  their  corresponding  normal  tissues  (Fig.  3A).  As 
many  collagens  and  collagen-remodeling  genes  are  frequently  upregulated  in  fibroblast 
activation  associated  with  inflammation  and  fibrosis  in  the  absence  of  cancer,  use  of  these  genes 
as  therapeutic  targets  in  cancer  could  be  problematic.  Analysis  of  expression  profile  datasets 
show  that  levels  of  COL11A1  mRNA  in  inflamed  colonic  tissue  from  inflammatory  bowel 
disease  and  fibrotic  lung  tissues  are  not  significantly  different  from  those  in  corresponding 
unaffected  colon  and  lung  and  that  COL11A1  levels  associated  with  colon  inflammation  and  lung 
fibrosis  are  minimal  and  markedly  different  from  those  associated  with  colon  and  lung  cancers 
(Fig.  3B).  In  contrast,  levels  of  ACTA2,  the  gene  encoding  the  prototypical  marker  of 
myofibroblast  differentiation,  aSMA,  is  expressed  at  similar  levels  in  cancers  and  inflamed  or 
fibrotic  tissues  (Fig.  3B). 
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Fig.  3.  Increased  expression  of  COL11A1  in  cancer  and  low  expression  in  normal  tissues, 
inflammation  and  fibrosis.  (A)  Comparison  of  COL11A1  mRNA  expression  in  normal  tissues 
and  corresponding  cancers.  Box  plots  for  two  different  platforms  (U133Plus2  and  U133A)  were 
generated  using  datasets  and  software  available  through  the  Gene  Expression  across  Normal  and 
Tumor  tissue  (GENT)  portal  (medical-genome.kribb.re.kr/GENT).  The  y  axis  shows  log2  mRNA 
levels.  Average  expression  levels  in  normal  tissues  and  cancer  tissues  are  indicated  by  vertical 
dotted  green  and  red  lines,  respectively.  (B)  COL11A1  and  ACTA2  mRNA  expression  in  normal, 
inflammatory  and  fibrotic  conditions  in  comparison  to  cancer.  The  graphs  were  generated  using 
the  public  R2  MegaSampler  software  (hgserverl.amc.nl/cgi-bin/r2/main.cgi)  for  the  processing 
and  normalization  of  individual  datasets  imported  from  the  Gene  Expression  Omnibus  (ul33p2, 
MAS5.0  platform).  The  number  of  samples  in  each  GSE  dataset  are  indicated  in  parentheses. 

3.3.4.  A  consistent  set  of  genes  is  co-expressed  with  COL11A1  across  different  cancers 

To  better  understand  the  biology  of  cancers  with  high  levels  of  COL11A1,  we  identified  genes 
that  most  closely  correlate  with  COL11A1  mRNA  expression  in  13  TCGA  datasets  representing 
different  cancer  types.  Spearman’s  rank  correlations  between  COL11A1  and  its  co-expressed 
genes  for  each  cancer  type  were  calculated.  The  genes  were  then  ranked  based  on  the  average 
correlation  of  each  gene  across  the  13  cancer  types.  The  top  195  correlated  genes  were  selected 
based  on  an  average  correlation  of  >  0.4  COL11A1  -correlated  genes  were  then  ranked  based  on 
the  average  of  the  absolute  correlation  values  (Table  1).  The  top  10%  most  highly  correlated 
genes  in  each  cancer  type  are  highlighted  in  pink  (Table  1).  Notably,  COL11A1- correlated  genes 
with  high  average  correlation  scores  also  tended  to  be  among  the  top  10%  highest  scored  genes 
in  each  cancer  type  (indicated  in  pink  in  Table  1). 
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1 92  CSGALNACT2 

193  BCAT1 

1 94  AXL 

195  CILP2 
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0.412 

0.39 

0.29 

0.55 

0.43 

0.53 

0.3 

0.33 

0.45 

0.52 

0.5 

0.34 

0.34 

0.37 

0.411 

0.54 

0.55 

0.69 

0.44 

0.5 

0.29 

0.22 

0.38 

0.64 

0.62 

0.19 

0.04 

0.23 

0.41 

0.62 

0.22 

0.73 

0.47 

0.66 

0.24 

0.42 

0.44 

0.75 

0.38 

0.28 

0.13 

-0 

0.41 

0.51 

0.47 

0.25 

0.35 

0.41 

0.36 

0.43 

0.5 

0.56 

0.55 

-0.1 

0.4 

0.63 

0.409 

0.59 

0.44 

0.24 

0.41 

0.32 

0.58 

0.44 

0.66 

0.18 

0.73 

-0.3 

0.35 

0.67 

0.409 

0.46 

0.27 

0.68 

0.38 

0.58 

0.44 

0.16 

0.48 

0.62 

0.49 

0.16 

0.23 

0.36 

0.408 

0.48 

0.55 

0.67 

0.24 

0.44 

0.57 

0.47 

0.63 

0.64 

0.48 

-0 

0.15 

0.01 

0.407 

0.37 

0.58 

0.18 

0.26 

0.48 

0.53 

0.3 

0.25 

0.63 

0.54 

0.24 

0.26 

0.66 

0.406 

0.6 

0.4 

0.51 

0.23 

0.7 

0.41 

0.25 

0.41 

0.5 

0.22 

0.09 

0.37 

0.58 

0.405 

0.28 

0.43 

0.71 

0.18 

0.25 

0.51 

0.48 

0.68 

0.28 

0.76 

-0.2 

0.45 

0.45 

0.405 

0.55 

0.46 

0.63 

0.55 

0.68 

-0.1 

0.16 

0.55 

0.65 

0.46 

0.2 

0.21 

0.24 

0.405 

0.47 

0.43 

0.64 

0.23 

0.42 

0.23 

0.33 

0.37 

0.52 

0.39 

0.37 

0.31 

0.55 

0.405 

0.52 

0.64 

0.76 

0.31 

0.42 

0.45 

0.24 

0.12 

0.31 

0.08 

0.19 

0.62 

0.58 

0.403 

0.55 

0.57 

0.66 

0.21 

0.38 

0.43 

0.14 

0.37 

0.54 

0.41 

0.17 

0.22 

0.57 

0.402 

0.28 

0.46 

0.26 

0.43 

0.46 

0.57 

0.3 

0.71 

0.5 

0.39 

0.19 

0.41 

0.25 

0.401 

Table  1.  COL11A1- correlated  genes  (Spearman’s  rank  correlation)  across  13  different  TCGA 
carcinoma  types.  Pink  denotes  the  top  10%  COL11A1  -correlated  genes  in  each  individual 
carcinoma  type. 
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In  contrast,  the  top  10%  COL11A1- anticorrelated  genes  were  not  conserved  across  these  cancer 
types  (data  not  shown).  Some  of  the  top  ranked  COL11A1- anticorrelated  genes  in  individual 
cancer  types  were  associated  with  normal  functions  of  these  organs  suggesting  that  they  may 
represent  normal  tissue  or  a  noninvasive  tumor  component.  For  example,  the  ovarian  cancer  top 
100  COL11A1- anticorrelated  genes  present  in  the  GSE12172  ovarian  cancer  dataset  were 
primarily  expressed  in  ovarian  tumors  of  low  malignant  potential  (Fig.  4). 


Serous  ovarian  tumors  GSE12172 
□  Low  malignant  potential  ovarian  tumor  (32) 
■  Invasive  ovarian  carcinoma  (58) 


Fig.  4.  Ovarian  cancer  COL11A1- 
anticorrelated  genes  are  enriched  in 
ovarian  tumors  of  low  malignant  potential. 
GeneSet  clustering  analysis  of  low 
malignant  potential  serous  ovarian  tumors 
(n=32)  and  invasive  serous  ovarian  tumors 
(n=58)  using  the  top  100  COL11A1- 
anticorrelated  genes  (Spearman’s 
correlation.  Table  S2)  from  the  ovarian 
cancer  TCGA  dataset.  The  clustering 
analysis  was  performed  using  R2 
(http://hgserver  1  .amc.nl/cgi- 
bin/r2/main.cgi). 


3.3.5.  Pan-cancer  COL11A1 -correlated  genes  are  induced  in  cancer-activated  fibroblasts 

Four  different  public  datasets  were  used  to  determine  the  cell  type  associated  with  expression  of 
the  pan-cancer  COL11A1- correlated  genes.  Several  of  the  195  pan-cancer  COL11A1- correlated 
genes  have  been  shown  to  play  a  role  in  EMT  [14].  Malignant  cells  undergoing  EMT  have  been 
proposed  as  one  possible  source  of  cancer  fibroblasts  [15].  To  determine  if  the  pan-cancer 
COL11A1- correlated  gene  set  is  preferentially  expressed  in  cancer  cells  undergoing  EMT  or  in 
host-derived  fibroblasts,  we  used  the  e-mtab-991  public  transcription  profile  dataset  of  primary 
patient-derived  colon  cancers  and  their  patient-derived  xenografts  (PDX)  in  nude  mice  [16]. 
Presumably,  in  PDX  samples,  fast-proliferating  human  cancer  cells  continued  to  grow  in  mice 
while  slow-proliferating  human  cancer  fibroblasts  were  lost  and  eventually  replaced  by  mouse 
fibroblasts,  which  can  be  distinguished  from  human  cells  by  species-specific  gene  probes  [16]. 
GeneSet  clustering  analysis  showed  that  most  of  the  pan-cancer  COL11A1- correlated  genes  had 
diminished  levels  in  PDX  samples  in  comparison  to  primary  cancers  (Fig.  5A),  indicating  that 
the  genes  are  enriched  in  the  fibroblasts  rather  than  in  the  cancer  cells.  The  preferential 
expression  of  the  COL11A1- correlated  genes  in  fibroblasts  was  further  confirmed  using  GeneSet 
clustering  analysis  of  the  GSE40595  dataset  in  which  ovarian  cancer  fibroblasts  and  epithelial 
cancer  cells  were  isolated  by  laser-capture  microdissection  [17].  The  pan-cancer  COL11A1- 
correlated  genes  were  preferentially  expressed  in  cancer  fibroblasts  (Fig.  5B).  In  addition  to 
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fibroblasts,  immune  cells  are  a  major  component  of  the  tissue  microenvironment.  To  exclude  the 
possibility  that  the  pan-cancer  COL11A1- correlated  gene  set  represents  immune  cells  in  the 
tumor  microenvironment,  we  used  the  expression  profile  of  230  mouse  hematopoietic  cell  types 
generated  by  the  Immunological  Genome  Project  (ImmGen)  compendium  [18].  In  addition  to 
immune  cells,  the  dataset  contains  expression  profiles  of  skin  fibroblasts  and  fibroblasts  residing 
in  the  thymus,  lymph  nodes,  and  spleen.  The  pan-cancer  COL11A1- correlated  gene  set  was 
highly  represented  in  fibroblasts  but  not  other  hematopoietic  cell  lineages  (Fig.  5C).  Cancer- 
activated  fibroblasts  have  a  different  expression  profile  than  normal  fibroblasts.  Moffitt  and 
colleagues  defined  a  23-gene  signature  of  ‘normal  stroma’  and  a  25-gene  signature  of  ‘activated 
stroma’  [19]  using  non-negative  matrix  factorization  for  virtual  microdissection  of  primary  and 
metastatic  pancreatic  ductal  cancer  samples  into  cell  subsets  with  prognostic  and  biologic 
relevance.  None  of  the  23  (0%)  ‘normal  stroma’  genes  in  contrast  to  18  of  25  (72%)  ‘activated 
stroma’  genes  were  present  in  the  COL11A1  -correlated  gene  set,  respectively  (Fig.  5D), 
suggesting  that  the  COL11A1- correlated  gene  set  represents  cancer-activated  fibroblasts.  To 
determine  whether  the  pan-cancer  COL11A1  -correlated  gene  set  is  associated  with  patient 
survival  in  the  -18,000  cases  of  liquid  and  solid  malignancies  in  the  PRECOG  dataset  [20],  we 
compared  survival  z-scores  for  the  195  pan-cancer  COL1 1A1  -correlated  genes  with  the  survival 
z-scores  for  all  genes  in  the  dataset.  This  analysis  showed  that  expression  of  the  pan-cancer 
COL1 1A1 -correlated  gene  set  is  significantly  associated  with  poor  survival  (Fig.  5E). 


A 


Colon  Cancere-mtab-991 
■  Prim  ary  turn  or  (40) 

□  PDX  (75) _ 
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B 


Ovarian  Cancer GSE40595 
□  Ovarian  cancer  epithelia  (32) 

■  Ovarian  cancer  fibroblasts  (31 ) 


D 


23  normal 

stroma  genes 
Moffitt,  2015 


195  pan-cancer 
COL11A1- 


correlated  genes 


25  activated 

stroma  genes 
Moffitt,  2015 


COL11A1- 

correlated 

genes 


Fig.  5.  The  pan-cancer  COL1 1A1  -correlated  gene  set  is  expressed  in  cancer-activated 
fibroblasts  and  associated  with  poor  patient  survival  in  multiple  cancer  types.  (A)  GeneSet 
expression  clustering  analysis  of  40  primary  colon  cancer  samples  and  75  patient-derived 
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xenograft  (PDX)  samples  in  the  e-mtab-991  dataset  using  the  pan-cancer  COL11A1- correlated 
genes.  (B)  GeneSet  expression  clustering  analysis  of  laser-microdissected  ovarian  cancer 
epithelial  cells  (32  samples)  and  fibroblasts  (31  samples)  in  high  grade  serous  ovarian  cancer  in 
the  GSE40595  dataset  using  the  pan-cancer  COL11A1- correlated  genes.  The  Euclidean  distance 
clustering  analysis  heatmaps  in  (A)  and  (B)  were  generated  using  the  public  R2  GeneSet 
Clustering  Analysis  tool  (hgserverl.amc.nl/cgi-bin/r2/main.cgi).  (C)  Expression  of  the  pan¬ 
cancer  COL11A1- correlated  genes  mapped  on  the  transcriptome  of  individual  murine 
hematopoietic  and  stromal  cell  types  in  the  ImmGene  project  (immgen.com).  The  plot  was 
generated  using  MyGeneset  tool  (rstats.immgen.org/MyGeneSet).  Black  dots  represent  mRNA 
levels  (y  axis)  of  186  pan-cancer  C0L11A1- correlated  genes  (9  genes  were  not  present  in  the 
database)  across  230  individual  cell  types  (X  axis)  grouped  into  10  main  groups.  (D)  Overlap  of 
the  195  pan-cancer  COL11A1- correlated  genes  with  23  ‘normal  stroma’  and  25  ‘activated 
stroma’  genes  defined  by  Moffitt  et  al.  (E)  Unweighted  meta  z-scores  of  191  COL11A1- 
correlated  genes  (4  genes  were  not  available  in  the  PRECOG  database)  compared  with  those  of 
all  genes  in  the  PRECOG  database.  The  plot  was  generated  using  GraphPad  Prism  software 
version  6.0.  Intergroup  differences  were  assessed  by  the  Student’s  t-test.  Mean  ±  SEM  of  pan¬ 
cancer  COL1 1A1  -correlated  genes  (0.8916  ±  0.1641  N=191);  Mean  ±  SEM  of  all  genes  (0.09918 
+  0.01416  N=23287);  ****p<0.0001. 
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4)  other  achievements 

■  What  opportunities  for  training  and  professional  development  has  the 
project  provided? 

Nothing  to  report. 

■  How  were  the  results  disseminated  to  communities  of  interest? 

Nothing  to  report. 


What  do  you  plan  to  do  during  the  next  reporting  period  to  accomplish  the 

goals? 

During  the  no-cost  extension  of  the  project,  we  will  validate  the  gene  signatures  in  patient 
samples  and  develop  a  preliminary  quantitative  assay  for  use  in  the  clinical  setting.  The 
development  of  a  reliable  test  for  the  identification  of  high  risk  patients  is  not  only  crucial  to 
improving  their  clinical  management  but  also  timely  because  of  the  emergence  of  personalized 
treatment  strategies  for  ovarian  cancer.  A  validated  gene  signature  to  identify  patients  with 
adverse  outcomes  has  the  potential  to  reduce  both  the  human  and  financial  costs  of  ineffective 
therapies  and  associated  toxicities. 
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4.  IMPACT: 

■  What  was  the  impact  on  the  development  of  the  principal  discipline(s)  of  the 
project? 

Future  efforts  to  specifically  target  activated  cancer  fibroblasts  in  ovarian  cancer  can  be 
improved  by  designing  novel  therapies  that  target  genes  such  as  COL11A1,  which  are 
specifically  expressed  in  activated  cancer  fibroblasts.  This  provides  an  opportunity  to  deliver 
targeted  therapies  directed  at  the  underlying  mechanism  of  the  poor  prognosis  signature.  This 
will  facilitate  more  individualized  treatment  decisions  and  improve  the  quality  of  care  for 
patients  with  ovarian  cancer.  In  contrast  to  aSMA  and  other  markers  frequently  used  in  mouse 
models,  COL11A1  is  not  expressed  in  mesenchymal  precursors  in  normal  organs  or  fibroblasts 
associated  with  non-cancerous  conditions,  such  as  inflammation  and  organ  fibrosis.  This  should 
increase  targeting  efficacy  and  reduce  off-target  toxicity. 

■  What  was  the  impact  on  other  disciplines? 

The  results  of  our  analysis  have  an  important  therapeutic  implication  for  solid  tumors  other  than 
ovarian.  The  strong  association  of  the  COL11A1 -coexpressed  gene  signature  with  poor  survival 
across  different  cancer  types  suggests  that  therapeutic  approaches  targeting  fibroblast  activation  in 
ovarian  cancer  should  be  effective  in  multiple  cancers. 

■  What  was  the  impact  on  technology  transfer? 

Nothing  to  report. 

■  What  was  the  impact  on  society  beyond  science  and  technology? 

Nothing  to  report. 


5.  CHANGES/PROBLEMS: 

■  Changes  in  approach  and  reasons  for  change.  No. 

■  Actual  or  anticipated  problems  or  delays  and  actions  or  plans  to  resolve 
them. 

Reasons  for  6-month  no-cost  extension  have  been  previously  reported  and  approved. 

■  Changes  that  had  a  significant  impact  on  expenditures.  No. 

■  Significant  changes  in  use  or  care  of  human  subjects,  vertebrate  animals, 
biohazards,  and/or  select  agents.  No. 

■  Significant  changes  in  use  or  care  of  human  subjects.  No. 

■  Significant  changes  in  use  or  care  of  vertebrate  animals.  No. 

■  Significant  changes  in  use  of  biohazards  and/or  select  agents.  No 

6.  PRODUCTS: 

Nothing  to  report. 
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■  Publications,  conference  papers,  and  presentations 

■  Journal  publications. 

Liu  Z,  Beach  JA,  Agadjanian  H,  Jia  D,  Aspuria  P-A,  Karlan  BY,  Orsulic  S.  Suboptimal 
cytoreduction  in  ovarian  carcinoma  is  associated  with  molecular  pathways  characteristic  of 
increased  stromal  activation.  Gynecologic  Oncology  2014;  139:394-400.  Lead  Article  with 
Editorial  Review  in  Gynecologic  Oncology;  Editorial  Review  in  Obstetrical  &  Gynecological 
Survey.  (Published,  acknowledged  grant  funding) 

Beach  JA,  Aspuria  P-J,  Cheon  D-J,  Agadjanian  H,  Walsh  CS,  Karlan  BY,  Orsulic  S. 
Sphingosine  kinase  1  is  required  for  TGF-P  mediated  fibroblast-to-myofibroblast  differentiation 
in  ovarian  cancer  Oncotarget  2015;  7:4167-4182.  (Published,  acknowledged  grant  funding) 

■  Books  or  other  non-periodical,  one-time  publications.  N/A 

■  Other  publications,  conference  papers,  and  presentations.  N/A 

■  Website(s)  or  other  Internet  site(s).  N/A 

■  Technologies  or  techniques.  N/A 

■  Inventions,  patent  applications,  and/or  licenses.  N/A 

■  Other  Products.  N/A 

7.  PARTICIPANTS  &  OTHER  COLLABORATING  ORGANIZATIONS 

■  What  individuals  have  worked  on  the  project? 


Name: 

Sandra  Orsulic 

Project  Role: 

PI 

Nearest  person  month  worked: 

1.2 

Contribution  to  Project: 

Dr.  Orsulic  oversaw  statistical  analyses  of  the  gene  signatures,  analyzed 
and  interpreted  the  data,  wrote  and  published  one  manuscript,  and 
prepared  another  one  for  publication. 

Name: 

Dongyu  Jia 

Project  Role: 

Postdoctoral  fellow 

Nearest  person  month  worked: 

3 

Contribution  to  Project: 

Dr.  Jia  organized  retrieval  of  pathology  samples,  isolated  RNA  and 
prepared  cDNA  from  samples,  optimized  quality  control  for  Nanostring 
assay,  and  assisted  in  the  writing  of  the  published  manuscript. 

Has  there  been  a  change  in  the  active  other  support  of  the  PD/PI(s)  or 
senior/key  personnel  since  the  last  reporting  period? 

Nothing  to  report. 
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■  What  other  organizations  were  involved  as  partners? 

Nothing  to  report. 

8  SPECIAL  REPORTING  REQUIREMENTS 

■  COLLABORATIVE  AWARDS:  N/A 

-  QUAD  CHARTS:  N/A 

9  APPENDICES:  N/A 
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